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ABSTRACT 

M star spectra, at wavelengths beyond 1.35 ^m, are dominated by water vapour yet 
terrestrial water vapour makes it notoriously dilRcult to make accurate measurement 
from ground-based observations. We have used the short wavelength spectrometer 
on the Infrared Space Observatory at four wavelength settings to cover the 2.5-3.0 
/xm region for a range of M stars. The observations show a good match with previ- 
ous ground-based observations and with synthetic spectra based on the Partridge & 
Schwenke (1997) line list though not with the SCAN (Jorgensen et al. 2001) line list. 
We used a least-squared minimisation technique to systematically find best fit pa- 
rameters for the sample of stars. The temperatures that we find indicate a relatively 
hot temperature scale for M dwarfs. We consider that this could be a consequence 
of problems with the Partridge & Schwenke linelist which leads to synthetic spectra 
predicting water bands which are too strong for a given temperature. Such problems 
need to be solved in the next generation of water vapour line lists which will extend 
the calculation of water vapour to higher energy levels with the good convergence 
necessary for reliable modelling of hot water vapour. Then water bands can assume 
their natural role as the primary tool for the spectroscopic analysis of M stars. 

Key words: stars: low-mass, brown dwarfs; stars: late-type; stars: abundances; stars: 
fundamental parameters; stars: atmospheres 



X 
S3 



1 INTRODUCTION 

More than two-thirds of stars within 10 parsecs are M dwarfs 
and it is very probable that this number density prevails 
throughout our Galaxy. Unless there is a sharp turn-down in 
the stellar mass function, they and even lower mass objects 
are a major component of the Galaxy's mass. The dominant 
red and infrared luminosity of the underlying stellar pop- 
ulation of galaxies is from M giants. The dominant source 
of opacity for late-type M dwarfs, giants and brown dwarfs 
is water vapour which easily forms in their relatively high 
pressure, low temperature atmospheres. Leaps in theoretical 
molecular quantum mechanics and computer hardware ca- 
pabilities mean that it is possible to perform ab initio calcu- 
lations to accurately predict the frequency and intensity for 
ro- vibrational transitions for water vapour. This means that 
it is no longer necessary to extrapolate laboratory measure- 
ments for water vapour to untestable temperature regimes 
of which are found in the atmospheres of M dwarfs. 

The preponderance of water vapour in the Earth's at- 
mosphere makes it very difficult to observe its spectral sig- 
nature in stars. At near-infrared wavelengths, where cool 
stars emit most of their flux, the strongest water vapour 



absorption band is centred around 2.65 microns where the 
atmosphere is opaque (Fig. 0). The advent of the Infrared 
Space Observatory has for the first time allowed observa- 
tions to be made at the peak of water vapour absorption 
in cool stars. Such data are not only impossible to obtain 
from terrestrial sites but also provides a vital overlap with 
ground-based data. For M dwarfs this is essential as the 
data reduction problems of decontamination of stellar and 
terrestrial water vapour are never far away. 



2 OBSERVATIONS 

The spectral region required for this program is inaccessi- 
ble from ground-based observatories and is well matched 
to ISO capabilities. SWS and ISO provide sufficient resolu- 
tion to resolve individual ro-vibrational water bands with 
enough sensitivity to observe the intrinsically faint and cool 
M dwarfs. The strategy was to observe a sample of M dwarfs 
which are bright enough to obtain high signal-to-noise spec- 
tra. The ISO programme for these observations is known as 
JONES_PROP32. 

We observed a range of M dwarfs - GJ 752B (M8V), 
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Figure 1. The water vapour absorption coefficient at 3000 K 
across the peak of the flux distribution for late-type M and brown 
dwarfs. The plot has been normalised to one around the maximum 
absorption of water vapour that can reliably be observed from the 
ground {'high opacity' case in Jones et al. 1994). 
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Figure 2. A comparison of datasets for GJ 406 from orbits 609 
and 717 of this programme and from UKIRT data from Jones et 
al. (1995). This object has the lowest signal-to-noise ratio in our 
analysis yet shows a reasonable agreement between spectra taken 
on different orbits and from the ground. 



Table 1. Properties and total integration times for the sample. 



Object 


Flux, Jy (2.76Atm) 


Spectral Type 


Time, s 


GJ 191 


4.13±0.09 


M2.5VI 


4530 


GJ 699 


6.66±0.14 


M4V 


1192 


GJ 406 


1.503±0.031 


M6V 


9839 


GJ 752B 


0.019±0.020 


M8V 


18216 


BS 8621 


509.7±4.1 


M4III 


1238 



GJ 406 (M6V), GJ 699 (M3.5V) and GJ 191 (M2VI) to- 
gether with the M4 giant BS 8621. The M dwarfs were cho- 
sen because they have been the subject of previous studies 
of M dwarfs (e.g. Jones et al. 1996). The observations were 
made in with the Short Wavelength Spectrometer (SWS) 
in its full resolution grating mode known as SWS06. Ob- 
servations were made using bands lA and IB covering the 
wavelength ranges 2.48-2.60, 2.60-2.75, 2.74-2.90 and 2.88- 
3.02 iim. These grating scans also gave simultaneous cov- 
erage from 15.99-16.25 and 16.03-16.21 /xm, however, the 
relatively low flux levels of the M dwarfs (<10 Jy) and the 
lower sensitivity of the instrument means that no useful data 
were recorded. 



3 DATA REDUCTION 

The observational data were processed within the Observers 
SWS Interactive Analysis Package from the standard pro- 
cessed data to Auto Analysis Result level through the stan- 
dard pipeline of Derive Auto Analysis Product. In addition 
to a few bad data points flagged by the auto analysis soft- 
ware ISO, SWS data also contains glitches and jumps not 
easily corrected for by the standard reduction tools (Heras 
1997). The philosophy adopted for the data presented here 
was to ignore all scans effected by glitches or jumps and 
to let individual bad points be taken care of by sigma clip- 
ping. The glitches and jumps do not tend to appear in the 
brighter objects though for the fainter objects, e.g., GJ 406, 



they affect around 10 percent of scans. The proc_band and 
div procedures were used to print out scans for each indi- 
vidual detector for each observation for all targets. These 
were visually inspected to check for glitches, and bad data 
points. The ISO Spectral Analysis Package (ISAP) was used 
to remove bad data points, check the data and to produce 
the final combined spectra by averaging across all detectors 
using the standard clip mean option. We found little sensi- 
tivity to the method used to combine the data. 

We also made observations of the archetypal late-type 
M dwarf GJ 752B. These data are tantalising however our 
conservative reduction philosophy means that little of this 
data survives. It would be possible to flag all glitches and 
jumps, however, they are usually accompanied by gradients 
affecting most of the data in a scan. We look forward to the 
advent of reduction tools that will be able to deal with flux 
levels at the 0.01 Jy level. 

The calibration files used for the reduction were those 
included in version 8.4 of ISAP package. Based on com- 
parisons between objects taken with the same configuration 
during different orbits the flux calibration varies by only a 
few percent. This is expected for band 1 SWS observations 
which have excellent calibration. The wavelength calibra- 
tion for the SWS instrument is measured to be within 1/8 
of a resolution element (Salama et al. 2000). The available 
ground-based data for our sample are at lower resolution 
than the ISO data, nonetheless, bear out the calibration 
of the SWS data. Fig. | shows the good agreement in the 
region of overlap between parts of the dataset taken on dif- 
ferent orbits and ground-based measurements (from Jones 
et al. 1995). 



4 AVAILABLE WATER LINELISTS 

The crucial opacity for M stars across the range of our 
observed spectra is water vapour (e.g.. Fig. 0). A useful 
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Figure 3. A spectral sequence for GJ 191 (M2VI), GJ 699 (M4V) and GJ 406 (M6V). 
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Figure 4. A comparison at optical spectral type M4 for the dwarf (V) GJ 699 and the giant (III) BS 8621. 



model atmosphere thus needs to include an appropriate wa- 
ter linelist. There are several different sources of water data 
and a few of them have been investigated here. A reliable 
water linelist needs: 1. a good (electronic) potential energy 
surface; 2. well converged nuclear motions (i.e vibration- 
rotation) calculations; 3. a reliable dipole surface. To date, 
it is not yet possible to a get a completely reliable ab initio 
potential surface, so all the linelists discussed in this pa- 
per used surfaces which have been adjusted to reproduce 
laboratory spectroscopic data for water. Fitting to lab data 
can cause problems in regimes where such data are unavail- 
able (Polyansky et al. 1997a). Conversely tests (Lynas-Gray, 
Miller & Tennyson 1995) have shown that ab initio dipole 



surfaces are much more reliable than the ones fitted to ex- 
perimental data. 

The water vapour data investigated in this paper are 
from Miller et al. (1995, known as the MT list), from Par- 
tridge & Schwenke (1997, known as the PS or AMES list), 
from Jorgensen et al. 2001 (known as SCAN), and from Viti 
(1997, known as the VT2 hnelist) which superseded the MT 
and VTPl line hsts (Viti, Tennyson & Polyansky 1997). MT 
is a relatively small list (10 million lines) which used, by to- 
day's standards, a rather inaccurate potential. VTPl energy 
surface is much more accurate than MT but it is still not 
complete enough for stellar models. VT2 is a much larger 
linelist ( > 300 million lines) which uses a potential surface 
reliable for higher vibrational states. However, due to seri- 
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Figure 5. Theoretical SEDs computed for 3100/5.5 model atmo- 
sphere taking into account continuum + atomic line and contin- 
uum + atomic lines -|- H2O opacities 
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Figure 7. Comparison between synthetic spectra at 3100 K using 
the SCAN and PS Une lists. 



Figure 6. GJ 406 compared with 3000 K synthetic spectra using 
the MT and PS linelists. The observed spectra and PS model are 
shown as solid thick and thin lines respectively and the MT model 
as a dotted line. The match of water bands between observations 
and model is considerably better for the PS model. 

ous computer format corruption problems this list could not 
be included in model atmospheres calculations. 



5 SPECTRAL ANALYSIS 

For this analysis we computed spectral energy distributions 
for a range of model atmosphere parameters: effective tem- 
peratures, Teff = 2700 — 3900 K, metallicities [Fe/H] = 

0.0 1.5 and gravities log g = 0.0 — 5.5. From here on 

we refer to these as model synthetic spectra. These were 
generated using the temperature structures from the model 
atmospheres of Hauschildt, AUard & Baron (1999) by the 
WITA6 program (Pavlenko 2000). In our computations we 
used the atomic line list from VALD (Kupka 1999) , the wa- 
ter vapour line lists of Partridge & Schwenke (1998) and 
Jorgensen (2001) and the CO line list of Goorvitch (1994). 
We also used synthetic spectra from the PHOENIX 



model atmosphere code, the so called NextGen version 5 
(Hauschildt et al. 1999). 

We adopt the Voigt funtion H{a, v) for the line shape 
of any line in our computations. Damping constants a — 
(72 + 74 + 76)/(4 X TT X Ai/o) are determined in atmospheres 
of cool dwarfs mainly by Van der Waals (pressure) broaden- 
ing {n—6) and (in upper layers) natural broadening {n—2). 
Stark broadening (n=4) of spectral lines can be neglected 
here due to low temperatures and electron densities. For 
atomic lines the line broadening constants were taken from 
VALD. Unfortunately, for molecular line broadening mod- 
els are rather uncertain. We assume that the mechanisms 
of broadening molecular and atomic lines are similar. Ra- 
diative broadening of molecular lines were computed in the 
frame of classical approach 72 — 0.22/ (Allen 1973). Due 
to the high pressures and low temperatures in the atmo- 
spheres of cool dwarfs Van der Waals broadening should 
dominate there (see Pavlenko 2001). The damping constants 
of pressure broadening of molecular lines 75 were computed 
following Unsold (1948). 

To provide an idea of the improvement since our previ- 
ous work (Jones et al. 1995), in Fig. ^we compare GJ 406 
with 3000 K synthetic spectra generated using the PS and 
MT lists. As with Jones et al. (1995) we found that the inten- 
sities of the water band strengths seem to be reasonably well 
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Figure 8. Number of energy levels as a function of energy for 
PS (dotted) and VT2 (continuous) linelist for J = 17. The PS 
linelist stops at 28000 cm~^. 

determined. However, we find that the PS line list predic- 
tion of bands with wavelength is a substantial improvement 
on the previously used MT line list. We also investigated 
using the SCAN line list (Jorgensen et al. 2001) and in Fig. 
1^ compare this to the PS list. As with our previous com- 
parisons to the SCAN list (Jones et al. 1996) we find the 
predicted shape of the water bands as well as the structure 
of the bands are rather different from other models and ob- 
servations. As expected from the scale of the calculation and 
previous comparisons with observations (e.g., AUard et al. 
2000) we find that the PS line list is the most useful water 
Une list available for our observed data. We first proceed to 
investigate the completeness of the PS line list and then use 
it to find best-fit parameters for our targets. 



5.1 Completeness of the PS line list? 

For model atmospheres with effective temperatures in the 
range 2000 — 4000 K, Boltzmann considerations suggest 
that transitions involving states with J = 20 — 30 are the 
dominant sources of opacity. We find two problems with the 
PS linelist: (i) first of all, it does not reach 30000 cm~^ 
for any J apart from J — 0. On average it gets to 28000 
cm~^. (ii) Secondly, even by truncating energy levels to 
28000 cm"\ VT2 stiU has more levels than PS has for high 
J, for example, at J = 17, VT2 gives 100 energy levels more 
than PS below 28000 cm~^. From the number of energy lev- 
els missing for each J we estimate that PS omits ~ 30 % 
transitions up to J = 10; ~ 50% up to J = 20; and ~ 60% 
up to J = 28. 

As a test, a new set of energy levels using the PS po- 
tential energy surface, the same atomic masses as PS and 
VT2 nuclear motion parameters, were computed. This new 
set of energy levels is known as UCLPS (University College 
London Partridge Schwenke list; see Polyansky et al. 1997b 
for further comparisons). This was done to test the varia- 
tional convergence for the same potential energy surface. An 
example of the energy levels comparison with VT2 and the 
new calculations is shown for J = 17 in Fig. ^. The calcula- 
tions performed with the PS potential energy surface are in 
close agreement with VT2 which suggests that the lack of 



states in the PS linelist is due to poor convergence and not 
to differences in the potential energy surface used. In fact, 
close examination of the parameters of the PS calculation 
strongly suggests that their decision not to increase the size 
of the Hamiltonian matrix beyond that used for J = 4 re- 
sulted in poorly converged calculations for higher states with 
high J. This decision undoubtedly saved them from some of 
the computational problems experienced in the computation 
of the VT2 linelist (Viti 1997). 

For low- J energy levels, PS's calculations give superb 
results, reproducing experiment with a much higher accu- 
racy than than VTPl or VT2. However, for higher rota- 
tional states, particularly those with J > 20, we find that 
a very high proportion of rotational states which ones ex- 
pects to be degenerate in fact show significant splittings in 
the PS linelist. This splitting is not shown in VTPl, VT2 
or in UCLPS. In particular for the levels with high Ka and 
with Kc odd (where Ka and Kc are the projections of J on 
the A and C principal axes of rotation of asymmetric top 
molecules), all lie below the ones with Kc even for levels 
with which they should be quasi-degenerate. Since the PS 
calculation truncates variational rotation-vibration calcula- 
tions with 7500 energy-selected basis functions independent 
of the rotational parity this means that the Kc odd calcu- 
lations will contain states of higher cut-off energy than the 
Kc even calculation. The variational principle means that 
the Kc odd states will be better converged and hence lower 
in energy. 

This causes two problems. Firstly, artificial splitting of 
lines means that it is difficult to use the list for line assign- 
ments (Polyansky et al. 1997a). The second is more subtle. 
An important consideration in radiative transport is how 
the line absorptions fill in gaps in the spectrum. Two tran- 
sitions which, to within their linewidth, are coincident will 
have less effect on the opacity than two separate transitions. 
Artificially doubling the number of lines for these J values 
is likely to cause the strength of water vapour bands at low 
resolution to be overestimated. 

It should be noted that the resolution of the observed 
spectra (see Fig. ^ is much too low to see the rotational 
structure of H2O shown in the theoretical spectrum (Fig. 
Although the calculated line positions of the PS list agree 
well enough with the experimental ones (tabulated, for ex- 
ample, in HITRAN96, Rothman et al. 1996), individual line 
strengths can differ greatly. Recently Schwenke & Partridge 
(2000) found problems with their analytical representation 
of the ab initio dipole moment data used for the computation 
of the PS linelist: this may have led to an overestimation of 
the intensities of weak bands with respect to experimental 
data. They improve their previous analytical representation 
by careful and accurate fitting. However their new dipole 
surface has yet to be used to produce an improved water 
linelist. 

5.2 Best-fit parameters 

Our synthetic spectrum models were calculated at high res- 
olution, typically 0.5 — 1 A. Before comparison with model 
spectra were made the observed data were corrected for their 
radial velocities: 244.0, 110.9 and 21.0 km/s respectively for 
GJ 191, GJ 699 and GJ 406 (Leggett 1992). The model were 
transformed to the resolution of the data and to wavelengths 
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Figure 9. An example of the dependence of S on /n and /g 
values of /n and /g are in A. 



GJ 406,3100/5.5, D„=0, f„=0.£ 




Figure 10. An example of the dependence of S on fs and /g, 
values of fs and /g are in A. 



in air. The resolution of ISO data are actually somewhat 
uncertain and have been the subject of a number of studies 
(e.g. Valentijn et al. 1996, van den Ancker et al. 1997, Lutz, 
Feuchtgruber & Morfill, 2000). The SWS06 mode is a grat- 
ing scan mode and in principle reproduces the full intrinsic 
resolution of the grating (so-called SWS02 mode). However, 
the SWS06 mode has poor sampling of a given wavelength 
element. For the purpose of the derivation of high-fidelity 
line profiles, this causes unwanted interactions between line 
profiles and the relative sensitivities for the 12 detectors of 
a band ('fiat-fielding in SWS-speak'), since all twelve de- 
tectors have to be combined for a good profile. This means 
that although the measured resolution of SWS02 in band 
lA (2.38 — 2.60 /im is 0.00Q7Mm and O.QQllAim in band 
IB (Lutz et al. 2000), the actual resolution of these SWS06 
spectra is likely to be somewhat less. 

In order to determine the best fit parameters for our tar- 
gets in a systematic fashion, we performed a least-squared 
minimisation on the spectra within a grid of synthetic spec- 
tra using the PS line list. As with the previous work we 
noticed that the form of the observed spectra along with its 
unresolved nature means that the line broadening employed 
by the models has a substantial influence on the fits. Al- 
though the SWS line profile is well determined to within 1 



per cent, the resolution of the data are not so secure (Lutz et 
al. 2000). To find the best fit model parameters we allowed 
the spectral normalisation, wavelength shift and resolution 
to vary. 

We compared observed fluxes H'^^ with computed val- 
ues U * Hl\7:. We let T/^bda = J {F'.-T * G{y) * dy, where 
the theoretical flux is F^^""'^ and G{y) is the instrumental 
profile modelled by a Gaussian. In our case G{y) may be 
wavelength dependent. To get the best fit we find the min- 
ima per point of the 3D function 

5(/n,/s,/g) = E(l-//=^"7//°'^=)^ 

We calculated this minimisation with each observed 
spectrum and the grid of synthetic spectra to determine a 
set of parameters fn (normalisation factor), fa (wavelength 
shift parameter), /g (resolution) for all the stars of our sam- 
ple. 

In general we found the dependence of S on /n and fs 
is much stronger than on the /g. For example see Fig. ^and 
px| for the case GJ 406. It should be noted that although 
the resolution across a grating setting is expected to be the 
same the resolution is expected to decrease slightly between 
the so called band lA (2.48 — 2.60 ^m) and band IB (2.60 
— 3.02 ^m) data. We account for this by also investigating 
the band lA and band IB data separately. Best fit values for 
/g were investigated for the observed spectra, across most of 
the grid of model atmospheres. They ranged from 7.0 (band 
lA) to 18.2 (band IB) A and always showed the expected 
resolution difference between band lA and IB data, e.g., ex- 
ample minimisation plots are shown for GJ 406 and a 3100 
K, log g — 5.5, solar metallicity model in Figs ^ and |l^. We 
generally found the largest and smallest values for /g when 
making comparisons with synthetic spectra with properties 
far from the best fit minimisation S values. For example in 
the case of GJ 406, the extreme best fit /g values of 8.0 
and 18.2 A are obtained with synthetic spectrum models of 
3300 K, [Fe/H]=-1.5 and 2700 K, [Fe/H]=0.0 respectively. 
Apart from giving us an independent check of the resolu- 
tion, this procedure served to investigate the sensitivity of 
minimisation S values to resolution /g. Although different 
/g's give significantly different values of S we found that 
the best fit temperature for a given /g was stable. We also 
found S values for observed and synthetic spectra smoothed 
by the same effective values. Although this gave the same 
best fit minimisation synthetic spectra, this conservative so- 
lution involves smoothing the observed data and led to a 
decrease in the sensitivity of the best fit solution. Given the 
numerical stability of our best fits to resolution and the lack 
of evidence for any changes in the SWS resolution with time 
we felt confident using the adopting the resolution values 
of 0.0012 and 0.0017 fim for band lA and band IB data 
as given in the SWS Observer's manual (de Graauw et al. 
1996). 

Our minimisation procedure also yields normalisation 
values fn for the different spectral regions. We find that the 
band lA data are fainter than the band IB data by 1-2% 
using our best fit synthetic spectra. Given that this is gen- 
erally better than the formal signal-to-noise of our spectra 
and that this result is somewhat model dependent we thus 
consider the flux calibration of the data to be robust and 
do not apply any new flux calibration to the data. The min- 
imisation procedure also yields values for the wavelength 
shift fs- We found a shifts of 3.7-4.3 A between observed 
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Figure 11. Minimisation S for different synthetic spectra com- 
pared to GJ 406 spectrum. 
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Figure 12. Tlie observed spectrum of GJ 406 compared to a 
3000 K, log g = 5.5, [Fe/H]=-0.5 synthetic spectrum. 



and synthetic spectra which corresponds to around 0.3 of 
our data resolution. This shift maybe due to the wavelength 
calibration of the data though at least part may arise from 
inaccuracies in the water vapour line list used to generate 
the synthetic spectra. 

Our results are also sensitive to the rather uncertain 
damping parameters adopted for molecular lines. To test this 
dependence two models at 2700 K, log g = 5.5, [Fe/H]=-0.5 
with different damping constants were compared. The liter- 
ature for atomic hues (e.g., Gurtovenko & Kostyk 1989) tells 
us that that a factor of two might be the maximum expected 
for this uncertainty due to the difference between classi- 
cal and quantum mechanical computations. Although both 
atomic and molecular lines are damped osciallations the os- 
cillations are greater in molecules so a larger factor maybe 
appropriate. Nonetheless it is also possible that molecular 
lines are less sensitive to damping because the structure of 
molecular levels maybe more flexible to changes of the sur- 
rounding electric fields: electrons of molecular orbitals are 
less bound in comparison with atoms. There may be some- 
thing like a "soft (quasi-)adiabatic response". We conducted 
a rather conservative experiment and compared models us- 
ing Unsold's 76 and with 5x7b and found a flux difference of 
around 5%. For our comparisons we consider this is an upper 
limit as in all cases we are investigating weaker molecular 
bands, due to higher effective temperatures and/or lower 
metallities. Furthermore, we are investigating the spectral 
differences in terms of synthetic spectra at relatively low 
resolution and find the expected uncertainty due to molec- 
ular damping constants to be much less than 1% and so not 
significant for our results. 

The following subsections investigate our fits to the each 
observed spectra in turn. The best fit models are then com- 
pared with a compilation of values from the literature in 
Section 5.6. 



5.3 GJ 406 

As discussed above and as illustrated in Figs ^ and ^ we 
carried a number of different tests on GJ 406 in order to 
find our preferred model fits. We tended to use GJ 406 as 
our primary observed spectrum because of its late spectral 
type and strong water bands. In addition to the tests dis- 
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Figure 13. Minimisation 5" for different synthetic spectra com- 
pared to GJ 699 spectrum. 



cussed above GJ 406 has data taken during four different 
orbits which allows us to analyse the orbits separately, thus 
checking any time dependence problems with the dataset. 
We found differences in the minimisation values obtained of 
around 7% for the 2.48 — 2.90 /im data and around 20% for 
the 2.88 — 3.02 ^m data though the best fit models were all 
within 100 K. 

Fig. |l^ illustrates the minimisation, S for various differ- 
ent synthetic spectra with the observed GJ 406 spectrum. 
For this plot we have split up the fits into those for data 
shortward and longward of 2.6 ^m. The improved minimi- 
sation for the higher resolution <2.6^tm data can be seen. As 
can the rather complex sensitivity of the minimisation with 
temperature and metallicity. The best fit value appears to 
be for [Fe/H] = -0.5 and 3000 K with the solar metallicity 
models suggesting 3100 K and [Fe/H] = -1.0 models sug- 
gesting 2900 K. An illustration of the quality of the fit can 
be seen in Fig. |l^ 

5.4 GJ 699 

GJ 699 (Barnard's star) is the largest proper motion object 
known outside our solar system and is thus suspected of 
being relatively old and having a low metallicity. In fact it 
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Figure 

3300 K, 



2.6 2.7 2.8 2.9 

Wavelength in microns 

14. The observed spectrum of GJ 699 compared to a 
log g = 5.5, [Fe/H]=-0.5 synthetic spectrum. 




Figure 

3700 K, 



2.7 2.8 2.9 

Wavelength in microns 

16. The observed spectrum of GJ 191 compared to a 
log g = 5.5, [Fe/H]=-1.0 synthetic spectrum. 
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Figure 15. Minimisation S for different synthetic spectra com- 
pared to GJ 191 spectrum. 
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Figure 17. The observed spectrum of BS 8621 compared to a 
2800 K, log g = 3.0, solar metallicity synthetic spectrum. 



tends to skew the proper motion properties and characteris- 
tics of the local neighbourhood so much that it is sometimes 
ignored from such studies. Nonetheless it is the most nearby 
and brightest known M4V dwarf and thus was included as 
target for this programme. Fig. ^ shows the minimisation 
plot for the observed spectrum. This suggests a best fit syn- 
thetic spectrum with 3300 K and [Fe/H]=-0.5. In Fig. 
we compare this fit with the observed spectrum. 

5.5 GJ 191 

GJ 191 has a high space velocity and displays CaH bands 
in its optical spectra and unlike GJ 699 is indisputably 
spectrally typed as a subdwarf. Our analysis oi S = 
/([Fe/If], Toff), e.g., Fig. ^ showed a best fit temperature 
of 3700 K and [Fe/H] values in the range -1.0 1.5. Com- 
putations with larger/smaller [Fe/H] provide systematically 
larger S. Unlike GJ 406 and GJ 699 we found that the match 
of observed and synthetic spectra was considerably better for 
data longward of 2.7 ^m. We suspected that this might be 
caused by the interference of CO Ai/ = 2 bands, so we inves- 
tigated higher (log iV(C) = -3.13) and lower (log N{C) = - 
3.68) carbon adundances than the scaled solar value (log 
N{G) = -3.48), however, we found no evidence that the dis- 



crepancies are caused by CO. Since the incompleteness of the 
water line list increases with increasing temperature (e.g.. 
Fig. ^ it is to be expected that the spectrum of GJ 191 is 
less well fit than the cooler dwarfs. AUard et al. (2000) reach 
a similar conclusion using the PS line list: 'The introduction 
of the PS-H20 opacities brings solid improvements of the 
near-infrared SED of late-type dwarfs but fails as the AH95 
models did to reproduce adequately the J-K colors of hotter 
stars'. Despite these problems a reasonable match between 
the observed and synthetic spectra can be seen in Fig. ^ 

5.6 BS 8621 

BS 8621 is a giant star with a well determined spectral type 
of M4ni. Unlike the cool dwarfs previously considered it is 
likely that the oxygen (and carbon) abundances in its atmo- 
sphere may be altered due to nuclear processing. Nonethe- 
less there are no indications in the literature that this well 
studied object has a composition different from solar so we 
start with the simplifying assumption that oxygen and car- 
bon abundances in the atmosphere of BS 8621 are are solar: 
log Ar(0)=-3.14 and log N(C)=-2,AS (Grevesse & Anders 
1989). 

At a given temperature we expect that water bands 
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Table 2. Best fit models for tlie sample and a collection of effective temperatures and metallicities from the 
literature. The value for BS 8621 is given without an error to reflect the uncertainty in this value. References 
are as follows (1) Krawchuk, Dawson & De Robertis (2000), (2) Gizis (1997), (3) Leggett (1992), (4) Cayrel 
de Strobel et al. (1985), (5) Bessell (1991), (6) Tinney et al. (1993), (7) Jones et al. (1994), (8) Leggett et 
al. (2000, apj, 535, 965), (9) Reid & Gilmorc (1984), (10) Veeder (1977) , (11) Pettersen (1980), (12) Mould 
(1976), (13) Basri et al (2000), (14) Kirkpatrick et al. (1993), (15) Tsuji et al. (1996), (16) Berriman et al. 
(1992), (17) Brett (1995), (18) van Belle et al. (2000). 

Values of T^ff (and [Fc/H]) from the literature (Reference given in caption) 

3466 -1.82(1), 3700 -1.5(2), -0.5 (4), 3550 -0.55 (12) 

2924 -1.35 (1), 3500 -0.5 (2), 3250(5), 3110 (6), 3095(7), 3100 -0.0 (8), 
3250 (9), 3110 (14), 3210 (15), 3130 (16), 3140 (17) 
2800(5), 2580(6), 2670(7), 2600 -0.0 (8), 2850 (9), 2800 (10), 2500 (11), 
2800 (13), 3000 (14), 2800 (15), 2565 (16), 2600 (17) 
3375-3475 (18) 



Object 


Best fit model 


GJ 191 


3700±100 K, 




[Fe/H]=-1.0±0.3 


GJ 699 


3300±100 K, 




[Fc/H]=-0.5±0.3 


GJ 406 


3000±100 K, 




[Fc/H]=-0.5±0.5 


BS 8621 


2800 K 



BS 8621 

2800/3.0/0, log N(C) = -3.48, log N(0) = -3.12 
log N(C) = -3.13 
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Figure 18. The observed spectrum of BS 8621 compared to 2800 
K, log g = 3.0 synthetic spectra with two different carbon abun- 
dances. The influence of an enhanced carbon abundance can be 
seen prominently on the CO bands between 2.49 and 2.52 /im. 



should weaken towards lower gravities as pressure broaden- 
ing lessens in importance. Although BS 8621 has a lower 
gravity and an early M spectral type, its spectrum has simi- 
lar spectral features although a somewhat different spectral 
shape. As for the M dwarfs its spectrum is clearly domi- 
nated by water vapour. However, the spectrum is markedly 
different on the left-hand side from the right-hand side. The 
longer wavelength side is well fit by a 2800K log g — 3.0 
model, see Fig. Though at wavelengths shorter than 
about 2.7 /im, water vapour absorption does not completely 
dominate the spectrum. It appears that transitions of CO 
Ai/—2 maybe important. In order to explain the < 2.7 ^im 
spectrum with CO features it is necessary to considerably 
increase the carbon abundance. This allows the CO features 
to be more clearly seen against the background of water 
vapour features. In Fig. a synthetic spectrum with a car- 
bon abundance of log (C) — 3.13 considerably improves the 
fit relative to a synthetic spectrum with a solar abundance 
pattern (log (C) = -3.48, log iV(0) = -3.12). From the lit- 
erature BS 8621 appears to be accepted as a close to solar 
metallicity standard M4III though Lazaro et al. (1991) find 
that BS 8621 has a reduced carbon abundance [C/H]=-0.5. 
Lazaro et al.'s result thus supports the expectation that the 



C/0 ratio in giants such as BS 8621 may be considerably 
modified by the CN-cyle. Such a drastic change in carbon 
abundance appears to have a relatively small effect on the 
synthetic spectra longward of 2.7 fim (Fig. p^ ) where CO 
bands are considerably weaker. In order to properly investi- 
gate this we need to use structural models computed for a 
variety of C/0 ratios. Such models are not currently avail- 
able to us and anyways are beyond the scope of the current 
paper where we wish to investigate the regime where water 
vapour is the dominant opacity, e.g. Fig. pi 



5.7 Discussion 

We have investigated a variety of model fits for a range of 
M stars. These are particularly promising for our cooler tar- 
gets. In Table ^ we present best fit parameters for our sam- 
ple though it is important to note that our study of the 
sensitivity of the spectra to the various model parameters 
is arguably more lasting than our determination of Tes, log 
g, and [Fe/H] values. We could have picked out these values 
much more quickly by eye though we preferred do this sys- 
tematically and thus be able to quantify the sensitivity of 
the spectra to different model parameters. 

In general the temperatures that we find are relatively 
high in comparison with the literature values which were 
obtained with a variety of different methods. An immediate 
concern is that our temperature for GJ 699 is higher than the 
3100 K obtained for the benchmark eclipsing binary system 
CM Dra composed of two M4.5 stars (Viti et al. 1997). The 
bolometric luminosity for GJ 699 is lower by around 20% 
so it would be expected to have a lower temperature than 
CM Dra. However, Viti et al. (2001) have recently demon- 
strated that CM Dra has a lower metallicity than GJ 699 
by around 0.5 dex. Based on Chabrier & Baraffe (1997) for 
a given luminosity this lower metallicity will cause a smaller 
radius by around 3% or more which may compensate for the 
luminosity difference allowing our temperature for GJ 699 
to be just acceptable. One possibility is that the effective 
temperature measured by spectra from 2.5 — 3.0 fim is 
not a good reflection of effective temperature of the spec- 
tral energy distribution. Since the opacity of water vapour 
dominates the infrared opacity budget and is particularly 
high in this region we would expect to see little contribution 
at r ~ 1 in this region from deeper hotter layers thus we 
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would actually expect to see cooler effective temperatures 
than expected but not hotter. A more plausible explanation 
for the high temperatures we find is the non-physical line 
splittings of the PS line list. These are likely to over predict 
the strength of water vapour transitions for a given temper- 
ature. This will lead to higher temperatures being fit to a 
given observed spectrum. 

We found very good fits for GJ 699 and GJ 406, how- 
ever, our fits for GJ 191 and BS 8621 are not so impressive. 
We consider that the relatively poorer fit for GJ 191 arises 
because of its higher temperature though we attribute our 
fitting problems with BS 8621 to the relative importance 
of CO in its atmosphere. For the case of GJ 191, our fit 
parameters are within the wide spread of parameters found 
by other authors. In contrast to the relatively high tempera- 
tures found for dwarfs our formal solution for the giant finds 
a significantly low temperature. We observe water vapour 
bands that are considerably stronger than are expected for 
its low gravity. Given that the temperatures for M4III giants 
are determined using interferometry results our 600 K lower 
temperature are notable! However, our fit to BS 8621 does 
not account properly for its apparent CO absorption bands. 
The evidence from Lazaro et al. (1991) is that BS 8621 is 
carbon deficient. A lower value of log A'^(C) is likely to dra- 
matically alter the chemical equilibrium, and as a result, 
CO and H2O opacities and the temperature structure of the 
model atmosphere. In particular the CO opacity will de- 
crease with a corresponding increase in H2O opacity. This 
increase in II2O opacity is likely to be important accross a 
wide range of effective temperatures. In the case of BS 8621 
this increase in H2O opacity will lead to an increase in the 
relative strength of water bands for a given temperature and 
thus cause us to fit BS 8621 to a higher temperature. Fur- 
thermore, as discussed above, we expect this spectral region 
to yield relatively low temperatures. This idea will soon be 
tested as van Belle et al. (2001) are investigating the differ- 
ent effective temperatures found when using effective radius 
measurements made at different wavelengths from V to K. 
Another part of the solution to resolve the temperature dis- 
crepancy is likely to be provided by the presence of a warm 
molecular sphere residing above the photosphere, the so- 
called MOLsphere (Tsuji 2000a & 2000b). However, given 
that we do see CO bands from 2.5-2.7 ^ra it is not possi- 
ble to properly test this scenario until spectral analysis has 
been done with appropriate structural models and a variety 
of C/O ratios. 



6 CONCLUSIONS 

Observations of water vapour show a good match with previ- 
ous ground-based observations and indicate that the PS line 
list predicts the positions and intensities of water vapour 
well enough to use for the determination of effective tem- 
peratures. The PS line list is a substantial improvement on 
the MT line list which has been widely used for the gener- 
ation of synthetic spectra of M dwarfs. The SCAN line list 
produced results substantially different from the observa- 
tions and other models. The effective temperatures we de- 
termined are reasonably consistent with other methods of 
temperature measurement. Although the complex spectral 
energy distributions of cool dwarfs caused by water vapour 



have traditionally hampered reliable temperature determi- 
nations, it now seems feasible that observations of water rich 
regions coupled with the next generation of water vapour 
line lists might become the method of choice for the tem- 
perature determination of cool dwarfs. Our work indicates 
that most sensitive best fits will be obtained when analysing 
high resolution data. 

However, a further conclusion of this work, in agree- 
ment with that of AUard et al. (2000), is that the presently 
available water linelists are still not good enough to generate 
accurate spectroscopic models of cool stars. These linelists 
have however been instrumental in improving the interpreta- 
tion of hot water, e.g. in sunspots (Polyansky et al. 1997b); 
as a result there are now significantly more experimental 
data available on water (Tennyson et al. 2001). These data 
are being used to greatly improve the effective potential en- 
ergy surfaces for water which, when combined with the im- 
proved dipole surface due to Schwenke & Partridge (2000), 
should make an excellent starting point for generating a new 
linelist. Experience, such as the tests performed in this pa- 
per, show that all aspects of such calculations need to be 
of very high quality if a satisfactory water opacity is to be 
obtained. 
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